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We study excitations of the spin-1 Bose gas at finite temperatures and in the presence of a not 
so strong magnetic field, or equivalently, when the gas sample is partially polarized. Motivated by 
the success of two-fluid hydrodynamics of scalar superfluids we develop a three-fluid hydrodynamic 
description to treat the low frequency and long wavelength excitations of the spin-1 Bose gas. We 
derive the coupled linear hydrodynamic equations of the three sounds and evaluate them numerically 
in a self-consistent mean field approximation valid for the dilute gas at the intermediate and critical 
temperature regions. In this latter region we identify the critical mode. 
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I. INTRODUCTION 

Multicomponent dilute gases are fascinating systems 
that can exhibit — at sufficiently low temperatures — 
phase transitions where the interplay between a super- 
fluid transition and the breaking of some other contin- 
uous symmetry can be observed. Maybe the simplest 
example is the Bose-Einstein condensation of the spin-1 
Bose gas [1, 2], where the system has a global SU(2) (spin 
rotation) symmetry as well as the usual U(l). For the 
weakly interacting gas Bose-Einstein condesation always 
infers the breaking of the SU(2) symmetry too, because 
the three component condensate wave function violates 
the spin rotation symmetry, even when spontaneous mag- 
netization is zero. Depending on the sign and magnitude 
of the spin dependent part of the interaction a magnetic 
phase transition can occur prior to or at the onset of 
Bose condensation [3-6]. Note, that for relatively large 
interactions it has been shown recently that a pairing 
transition occurs instead of Bose condensation [7]. 

In the Bose condensed phase the global U(l) symme- 
try is always broken and correspondingly a Goldstone 
mode has to exist among the elementary excitations of 
the system. Though the spin rotation SU(2) symmetry is 
also broken by the condensate wave function two possi- 
blities still remain. (1) A U(l) symmetry — the rotation 
around the 2:-axis — remains a symmetry, at least with 
the combination of the global U(l) gauge symmetry. Or, 
(2) even this kind of z-axis rotation symmetry is broken 
and the ground state has no remaining symmetry at all. 
In this latter case another Goldstone mode must exist. 
This kind of double symmetry breaking was illustrated 
in the Random Phase Approximation (RPA) for the po- 
lar spin-1 gas earlier, and the Goldstone modes have been 
identified amongst the quasiparticle excitations [8]. 

However, identifying the critical mode and studying its 
properties can be achieved in the hydrodynamic approxi- 
mation which can be contimied by scale invariance to the 
critical region. Such a hydrodynamic approximation, the 
celebrated 2-fluid hydridynamics of superfluids, is the ba- 
sis of understanding the critical properties of liquid ''He 



even around the A-transition and also serve for under- 
standing the behavior of dilute, scalar Bose gases when 
interparticle collisions take place sufficiently frequently 
[9-11]. In this paper we will follow the route originally 
developed for supcrfluid ^He [12], and derive a three- 
fluid model to describe the low energy excitations of the 
system. We will evaluate the thermodynamic quantiti- 
ties, entering to the hydrodynamic equations, in a self- 
consistent mean-field approximation that treats the non 
condensed atoms as an ideal Bose gas above the spinor 
BEG. 

The rest of the paper is organized as follows. In Sec. 
II. the model is introduced with the help of the grand 
canonical Hamiltonian and its symmetry properties are 
discussed. In Sec. III. the two possible phases and the 
phase diagram of the system is given in the magnetic 
field-temperature plane. Three fluid hydrodynamics is 
developed for the partially polarized phase in Sec. IV, 
with some details of the derivation moved to the; Ap- 
pendix. The thermodynamic quantities are evaluated in 
the Bogoliubov-Hartree approximation and are presented 
together with the results for the hydrodynamic excita- 
tions also in Sec. IV. Finally we summarize the results 
in Sec. V. 



II. HAMILTONIAN OF THE SPIN-1 BOSE GAS 



We consider a system of homogeneous, weakly inter- 
acting, dilute, spin-1 Bose gas at ultralow temperatures 
and in a homogeneous magnetic field. The interparticle 
interaction is modelled by s-wave scattering, i.e. we ne- 
glect the relatively weak dipolar interaction of the gas. 
The grand-canonical Hamiltonian of the system takes the 
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following form: 

"H = ^ [(e^: - iJ.)Srs - gfJ-BB (F^)^^ al{k)as{k) 
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where ar{k) and a\{k) are the annihilation and cre- 
ation operators of plane wave states with momentum 
k and spin projection r. The spin index r refers to 
the eigenvalue of the z-component of the spin opera- 
tor and can take values from -t-,0,—. Correspondingly 
Fz = diag(l, 0, —1) is a 3x3 diagonal matrix. In Eq. (1) 
= fi^/c^/(2M) refers to the kinetic energy of an atom, 
fjL to the chemical potential, g to the gyromagnetic ratio, 
//B to the Bohr magneton, B to the modulus of the ho- 
mogeneous magnetic field. V is the volume of the system 
and V^g'^ the amplitude of the two particle interaction, 
given for spin-1 bosons by [1, 2, 13]: 



v. 



(2) 



with Cn = 47rfi2(ao + 2a2)/(3Af) , and Cs = 4:ntf{a2 - 
ao)/(3M). The parameters ag and 02 are the scattering 
lengths in the total hyperfine spin channel zero and two, 
respectively. The constant c„ > 0, while can both be 
positive or negative, depending on the relative values of 
oo and 02- If Cs < the interaction tends to align the 
spins, while for > a zero net spin is energetically 
favorable (in the absence of a magnetic field). For this 
reason systems with < are referred as ferromagnetic 
systems, while those with Cg > are called polar gases 
[1]. For example the ultracold gas of ^''Rb atoms in the 
f ~ 1 hyperfine state is ferromagnetic [14], while the 
^■^Na gas (also in the f — I hyperfine state) is polar [15]. 

The Hamiltonian (1) is invariant under the global U(l) 
symmetry transformation, which leads to particle num- 
ber conservation. In the case, when the external mag- 
netic field is zero, the Hamiltonian is also invariant un- 
der the global spin rotation SU(2), which leads to the 
conservation of the total magnetization (spin) of the sys- 
tem. When a magnetic field is present (breaking explic- 
itly the rotational invariance), only another U(l) sym- 
metry remains, namely the rotation around the z-axis. 
These two U(l) symmetries can be described either by 



a± — > a-|-e"^± together with ao 
equivalently by 
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or 



It is easy to see that 
(j) = + tf-) corresponding to the symmetry gener- 

ated by the particle number, and — — </?-) corre- 

sponding to the symmetry responsible for the conserva- 
tion of the 2-component of the total spin. The constraint 
due to the conservation of the total magnetization can be 
resolved similarly as that of the conservation of particle 
number, i.e. with the introduction of a Lagrange multi- 
plier in the Hamiltonian (1). This multiplier shows up 
in the same way as the magnetic field does, therefore an 




FIG. 1. (Color online) The phase diagram of the polar spin-1 
Bose gas in the Bogoliubov-Hartree approximation. -Bc(O) is 
the critical magnetic field at zero temperature, above which 
the condensate is fully polarized. Here we have chosen the 
coupling constants as c„ = 0.2 x MknTc/p and Cs — 0.4 x c„. 



effective magnetic field can be introduced as a sum of the 
external magnetic field plus the Lagrange multiplier. In 
the following B will mean this kind of effective magnetic 
field even when the physical magnetic field is zero. In 
the following we assume that either such is the case, or 
the external magnetic field is so small that the quadratic 
Zeeman effect can be neglected. 



III. PHASES OF THE SPIN-1 POLAR BOSE 
GAS 

The weakly interacting spin-1 Bose gas can have two 
types of Bose-Einstein condensed phases. One is the fully 
polarized Bose-condensed phase that is realized by the 
thermodynamic ground state of a gas with Cg < (fer- 
romagnetic coupling), or even when Cs > but the gas 
sample is originally polarized (or equivalently when the 
(effective) magnetic field is sufficiently large) . This phase 
will be referred to as the PI phase, in analogy to the 
similar phase of '^He. In this phase the condensate wave 
function is single component, and correspondingly only 
one of the two global symmetries are broken, namely the 
symmetry parametrized by the angle. In the case of 
a polar gas and starting from the PI phase, when the ef- 
fective magnetic field is lowered, the system undergoes 
a "second" Bose-Einstein condensation. The conden- 
sate wave function is no longer single component and the 
other remaining global symmetry becomes broken too [8] . 
This phase will be called as the P2 phase and one can as- 
sume that the condensate wave function has two nonzero 
components. We illustrate the phase diagram of polar 
spin-1 gas in Fig. 1 in the self-consistent Bogoliubov- 
Hartree approximation, described later. Note, that the 
point (T = Tc,B = 0) is a special one, where two crit- 
ical lines cross each other, and a renormalization group 
analaysis around 4 — e dimensions to leading order in e 



3 



show that there is no stable fixed point [16] . 



IV. THREE FLUID HYDRODYNAMICS IN THE 
P2 PHASE 

A. Equations of motion 

The system has two conserved global charges associ- 
ated with the symmetry of the internal states, namely 
the total number of particles and the ^-component of the 
total magnetization, therefore two continuity equations 
arise 



dtp + Vffp = 0, 
dt(T + fga = 0, 



(3a) 
(3b) 



where p is the mass dcnisty and a is the spin denisty 
(multiplied with the mass of the particles M). Corre- 
spondingly, Qp is the total current (momentum) density 
and (jcr is the spin current density. 

In the P2 phase the Bose-Einstein condensates has two 
components, one is in the -|- and the other is in the 
— spin components. Directly following from the corre- 
sponding Gross-Pitaevskii equations and from Galilean 
invariance we have two independent superfluid compo- 
nents with their velocities obeying 



MdtVs,+ + V{n + 6) = 0, 
MdtVs,- + V(m - 6) = 0. 



(4a) 
(4b) 



The effect of fluctuations in the chemical potential /i, and 
in the magnetic energy b = QHbB is nicely exhibited in 
an alternative representation. Instead of Vs,+ and Vg,- 
one can introduce their linear combinations 



Vn = 



Va = 



(5a) 
(5b) 



Clearly Vp is the superfluid velocity for the density cur- 
rent and v„ is the superfluid velocity for the spin current. 
They obey the following equations of motion 



with s being the entropy density of the gas. 

To close the set of hydrodynamic equations one has to 
relate the complete density- and spin currents to the 3 
components of the superfluid, i.e. to the normal compo- 
nent and to the two distinct superfluid ones. The only 
linear relation consistent both with the Maxwell relations 
detailed in the Appendix and with the zero temperature 
results one has to choose them to be 



Qp = PnVn + PsVp + (TsVcr, 
9(7 = O-nVn + CTsVp + PsVa, 



(9a) 
(9b) 



with p = + Ps and a = an + os- One can also split Ps 
and (Tg according to 



Ps = Ps,+ + Ps 
CTs = Ps,+ - Ps 



(10a) 
(10b) 



where ps,+ and ps,- arc the superfluid densities in spin 
components -|-, and — , respectively. 



B. The wave equations 

With the help of the equation of continuity for the 
density (3a) and with the Euler equation (7) one can 
derive the first of the three coupled wave equations: 



dfp-V^p = 0. 



(11) 



The second equation follows from the continuity equa- 
tion for the spin density (3b) by combining it with the 
Euler equation (7) and the equations of motion for the 
superfluid velocities (6) and the Gibbs-Duham relation. 
It reads as 



P Pn 



(^{(TnPs - Pn(^s) - P{(^n(^s - PnPs) ^ q 

PPnM ' ^ ^ 

where we have introduced the s = s/p, the entropy per 
unit mass. The final equation is obtained from the en- 
tropy continuity equation (8) in combination with the 
Euler equation (7) and the equations of motion of the 
superfluid velocities (6): 



MdtVp + Vp = 0, 
Mdtv^ -h V6 = 0. 



(6a) 
(6b) 



Momentum conservation yields the well known Euler 
equation: 



dt9p + Vp = 0, 



(7) 



with p being the pressure field. Another equation of mo- 
tion is provided by the entropy transport of the normal 
component of the fluid: 



dtS + V{sVn) = 0, 



(8) 



dfS — s 



^2 Ps ^2 J, _ - 



-cr ps - pcr s^2 

Pn PPnM 



V^b = 0. 



(13) 



Close to equilibrium the fluctuation of the three den- 
sities can be expressed with the help of the fluctuations 

of the intensive parameters 



dp 

da 



9^P={ Z ( ^ ]0!t+{ ^ ]d^b, (14a) 



dp 



dtp- 



drs={ 



dp 
df 

d 



db 

da 

m> 

ds 



dfb, (14b) 



(14c) 



4 



Combining Eqs. (11), (12), (13), and (14) one gets a 
closed set of wave equations for the quantities p, T and 
b describing density waves, temperature waves and mag- 
netization waves. When looking for the plane-wave solu- 
tions one uses 





i{kf—ck t) 



(15) 



tained from the eigenvector problem: 




(16) 



and the plane wave amplitudes and dispersions are ob- where 



M 



(i). 



(If)' 



(if)' 



'^{lTnPs-Pn(ra)-p{<T„tT3-PnPa) 



db 



pp„M 
' PPnM 



(17) 



The dispersion follows from det M = which is a cubic 
equation for and gives the velocities of the first sound, 
the quadrupolar mode and the second sound. 



where A = ■\/27rfi^/ (Mfc^T) is the de Broglie wavelength, 
and C(s) is the Riemann-zeta function. The total densi- 
ties then 



C. Evaluation of the thermodynamic quantities 

In experiments made with ^^Na atoms the typical ratio 
of the coupling constants is e = Cs/c„ « 3 x 10~^ <C 1. 
In this limit the cubic equation for the speed of sound 
resulting from the vanishing of the determinant of Eq. 
(17) factorizes into a quadratic equation: the first and 
second sound modes, and to a separate linear equation 
with solution in the order of e: the quadrupolar spin- 
wave mode. This kind of separation is a generic feature 
of the dilute Bose gas and results from the behavior of 
the thermodynamic derivatives appearing in Eq. (17). 
We arc going to illustrate this feature in the Bogoliubov- 
Hartree approximation but keep the value of e general in 
order to study the hybridization of the different sound 
waves. 

The advantage of the Bogoliubov-Hartree approxima- 
tion is that it is an extension of the Bogoliubov approx- 
imation to finintc temperatiires with accounting inter- 
action effects with noncondensed atoms in a way which 
leads to a conserving and gapless approximation together 
with a continous phase transition [8]. 

In this approximation the chemical potential and the 
magnetic field is expressed as 



b = c 



'M' 
a 



(18a) 
(18b) 



The normal density of spin component r reads as 

MC(3/2) 



3MC(3/2) ^ 3/2 



Pn = Pn,+ + Pnfl + Pn - = ^3' ' = Pt^'^ (20a) 

Ps = Ps,+ + Ps,- = P(l - t^^^) , (20b) 

0-„ = pn,+ - Pn - = 0, (20c) 

Mb 

(^s = Ps,+ - Ps- = = , (20d) 



with t = T/Tc is the reduced temperature. Tc is the tem- 
perature of Bose- Einstein condensation of the spin-1 Bose 
gas in zero magnetic field (or with zero magnetization). 



27rfi2 / p 
ksM \2,C{2,/2)M 



2/3 



(21) 



The grand canonical thermodynamic potential in this 
approximation is 



$(T,y,M) = -^3C(5/2)-^(c„p2 + c,a2). (22) 



2M2' 



A3 



(19) 

The pressure p = —'^ /V , and its thermodynamic deriva- 
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FIG. 2. (Color online) The speeds of sound versus temperature at B = (a). The relative values of the three components of 
the eigenvectors corresponding to the given hydrodynamic mode: the first sound (a), the second sound (b) and the third sound 
(c). The coupling constants are chosen to be the same as in Fig. 1. 
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D. Results and discussion 



dp 
df 

dp 

dp 

dp 
db 

da\ 

df) 

da 

dp 

da 

db 



T,b 



p,b 



T,b 



a 

M' 



= 0, 



= 0, 



M 

Cs ' 



(23a) 
(23b) 
(23c) 
(23d) 
(23e) 
(23f) 



The specific heat under constant volume is 



cv = T 



(24) 



The three solutions of det M = provides the three 
speeds of sound of the hydrodynamic modes, which we 
will call here as first, second and third sound. To iden- 
tify which speed corresoponds to which physical exci- 
tation mode one has to look for the eigenvectors, i.e. 
the solution of Eq. (16), corresponding to the specific 
speed of sound substituted into the matrix M. In gen- 
eral the hydrodynamic modes hybridize, and their roles 
exchange. In Fig. 2 (a) we plot the three speeds of sound 
at the whole temperature range for zero magnetization 
(or equivalently at zero magnetic field). At this specific 
choice of zero magnetization the phase transition hap- 
pens at Tc and the P2 phase goes directly to the normal 
phase (compare with Fig. 1 at i? = 0). This transition is 
characterized by the simoltaneous vanishing of the con- 
densates both in r = -f and in r = — . Correspondingly 
two Goldstone modes are present, which are the second 
and third sounds. The speed of the first sound stays 
finite at this transition. On subfigures Fig. 2 (b)-(d) 
we demonstrate the amplitudes of the fluctuations of the 
sound modes. It can be seen that the first sound is ba- 
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FIG. 3. (Color online) The speeds of sound versus temperature at _B = 0.6Bc(0) (a). The relative values of the three components 
of the eigenvectors corresponding to the given hydrodynamic mode: the first sound (a), the second sound (b) and the third 
sound (c). The coupling constants are chosen to be the same as in Fig. 1. 



sicly a pressure wave when we are far away from the point 
of the hybridization of the first and second sounds. The 
second sound is mainly temperature wave after the point 
of crossing with the third sound and not too close to the 
hybridization with the first sound. The third sound is 
a typical magnetic field wave after the crossing with the 
second sound. It is also very nicely demonstrated that 
the nature of the second and third sound changes at the 
crossing point. 

On Fig. 3 the same quantities are plotted but for finite 
magnetic field (or finite magnetization) . It is visible from 
the figure, that it is the quadrupolar mode that goes soft 
at the P1-P2 transition point, thus it plays the role of 
the Goldstone mode. The second sound stays finite, since 
it only softens when the BEG gives way to the normal 
phase, that is only at the Pl-normal transition point, 
which is not discussed in this paper. The first sound 
remains unmodified by the change of the magnetic field, 
since it is basically a pressure wave. The second and third 
sounds are hybridizing in this case and the crossing is 
resolved to an avoided one in the presence of the magnetic 
field, the explicit breaking of the rotational symmetry. 

In the limit when e = Cg/cn <C 1 the pressure and 



temperature waves separate from the quadrupolar spin 
wave. This limit can be treated by looking for solutions 
of det M = in the order of e. To leading order this gives 
the following solution: 



M 



1 - <3/2 



(25a) 



which is the dispersion of the quadrupolar spin wave. For 
the other solutions one looks for c ^ 0(1), and to leading 
order gets 



ksTc 
M 



7c+ 



2 sM\ 2 , 2 fkBTr 



M 



1 



0, 

(25b) 

which is quadratic in and its solution gives the first 
and second sounds which formally agrees with that of 
the scalar Bose gas [10]. Note, that in the t limit 
we get back the solutions — ksTc-jc/M and = 0. 
Thus at low temperature the lower of the two solutions 
crosses the quadrupolar mode of Eq. (25a) , which stays 
finite at zero temperature. We have introduced the fol- 
lowing dimensionless parameters 7^ = Cnp/ {MkBT^.), 
u = Mb/ (csp) = a I p. 
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We have evaluated the thermodynamic derivatives in 
the Bogohubov-Hartrce approximation, which gives even 
quantitatively good interpolation between zero tempera- 
ture and the critical point as far as equilibrium proper- 
ties are concerned. As a matter of fact at intermediate 
temperatures also for the scalar gas we have evaluated 
the two fluid hydrodynamic equations in the Bogoliubov- 
Hartree aproximation and found the speeds of sounds 
very close to (few %) the Popov approximation discussed 
in detail in Ref. [10] . In the very low temperature region, 
however, where fc^T < /x, the thermodynamical quanti- 
ties are dominated by phonon excitations and show a 
qualtitively different behavior. The big difference is that 
the speed of the temperature wave does not go to zero but 
rather to c| = cf /3. Therefore neither of the three sounds 
vanishes at zero temperature. In the limit of e <C I the 
quadrupolar mode is always situated below the tempera- 
ture wave; thus the third sound is always a quadrupolar 
spin wave, while the second sound is a temperature wave. 

One can assume that the susceptibility given in 
Eq. (23f) remains singular when — >■ beyond 
the Bogoliubov-Hartree approximation. Then the term 
c^{da/dh) can remain finite if is of 0{cs), and in the 
matrix (17) all terms containing can be neglected in 
this limit except the one multiplied by the susceptibil- 
ity (first term of M2_3). As a result one arrives at the 
velocity of the quadrupolar spin mode in this limit as 

where Eqs. (10) have been used. One can easily convince 
oneself that approximating the supcrfluid densities by the 
corresponding condensate densities (26) agrees with the 
expression (25a). However, Eq. (26) is valid to any orders 
in c„ . By the same token one can write the generalization 
of Eq. (25b) as the solution of the equation det M2,3 = 0, 
where M2,3 denotes the submatirx obtained from (17) by 
removing the second row and the third column. 

Concerning the dynamical critical properties, the 
P2— )-Pl transition on the T,B > plane belongs to 
the same universality class as the A-transition of liq- 
uid Helium [17, 18]. Consequently the quadrupolar spin 
wave suffers a damping Dk^, where D oc A~^/^ when 
approaching the transition line between the P2 and PI 
phases along the thermodynamic path on which A mea- 
sures the distance from criticality. 



V. SUMMARY 

In this paper we have developed a three fulid hydro- 
dynamical approach for studying the low energy excita- 
tions of the spin-1 Bose gas in the partially polarized 
Bose- Einstein condensed phase (P2). Such a low energy 
and low momentum approximation of the excitation spec- 
trum is considered to be valid, when interparticle scat- 
tering events happen sufficiently frequently in order to 



have large enough regions in the gas sample where ther- 
modynamic equilibrium can be assumed. This kind of 
approximation is a generalization of the two fluid hydro- 
dynamics, which was applied to scalar Bose-Einstein con- 
densates successfuly [10, 19] and is in sharp contrast to 
quantum hydrodynamic approximations, already stiidied 
for spinor Bose-Einstein condensates too [20, 21]. These 
latter aproaches are nonlinear equations, equivalent to 
the Gross-Pitacvskii equation and are valid only at (or 
very close to) zero temperature and where quasiparticle 
interaction can be assumed to occur very rarely. 

In the partially polarized Bose-Einstein condensed 
phase of polar spin-1 Bose gases (P2 phase) the fluid 
mass density can be decomposed to three components: a 
normal mass density and two supcrfluid densities. The 
relative motion of these three fluids determine the low 
energy excitations of the system. We have calculated the 
three sound velocities governing these fluctuations. We 
have found, that the speed of sound of the magnetization 
wave is always below that of the temperature and pre- 
sussure waves for sufficiently weak spin-spin interaction 
(cg <^ c„). However, if Cg is in the same order as c„, 
successive hybrydizations occur, and the nature of the 
excitation branches exchange their character. One can 
hope that with the help of magnetic and optical Fesh- 
bach resonances such hybridizations can be achieved ex- 
perimentally, and the three fluid nature of polar spin-1 
Bose condensates and its implications on the phase tran- 
sition can be examined. 
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Appendix A: Thermodynamic quantities 

The internal energy per unit volume of the supcrfluid 
spin-1 system in the P2 phase should look like 

de = MT ds + fi dp + b d<7 + Vn dgp + j+ dv+ + j- dv- . 

(Al) 

The quantities v+ = {n/M)Vip+, w_ = {n/M)V(p- are 
the supcrfluid velocities of the -|- and — spin components, 
respectively, while j+ and j_ are their conjugate flelds. 
In the PI phase the last term is ommited from Eq. (Al). 
In an alternative treatment (j> and 9 can be introduced 
instead of 0-1- and (;6_ . It has an advantage of sonic phys- 
ical insight. In this formulation the internal energy can 
be given by 

de = MTds+iidp+bda+Vndgp+jpdvp+jadva- (A2) 
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To obtain further thermodynamic equations one has 
to consider the transformation properties regarding an 
infinitesimal Gaheleian transformation, dw. The nonin- 
variant quantities transform as follows: 



de = gp dw, 
dgp = pdw, 
dv± = dw. 



(A3a) 
(A3b) 
(A3c) 



With their help, the momentum density can be identified 
as: 



9p = pv„ + j+ + j- 



(A4) 



In the alternative representation of the superfluid ve- 
locities, according to Eqs. (5) it is clear, that one of 
the superfluid velocities introduced in Eq. (A2) trans- 
form like a velocity: dvp = dw, however the other one is 
Galilean invariant: dv„ = 0. Therefore the momentum 
density is as follows: 



Qp = PVn+jp- 



(A5) 



Close to thermodynamic equilibrium the currents can be 



expanded to linear order in the velocities: 

dgp ^ 



9p 



dv„ 



dVp'''' 



dvo 



dVcr 



dvn ovp dv^ 



(A6a) 
(A6b) 
(A6c) 

Equation (A2) yields the following Maxwell relations: 

(A7a) 
(A7b) 
(A7c) 

The nine coefficients of the linear expansion (A6) are 
not independent. With the help of the Maxwell realtions 
(AT) and Eq. (A5) 



dgp 


_ djp 


dvp 


dVn 


dgp 


_ _dja 


dVa 


dv„ 


djp 


_ dja 


dVa 


dvp 



jp = A{Vp-Vn) + BVa, 
Ja = B{Vp - Vn) + CVa 



(A8a) 
(A8b) 



hold for linear order. Choosing A,B,C to be consistent 
also with the zero temperature linearized hydrodynamic 
equations gives 



A 
B 
C 



djp 


- 


dvp 


dVn 


djp 


_ dja _ 


dv„ 


dVp 


dja 




dVa- 


= Ps- 



dgp dgp 
dvn 



dVn 



= Ps, 



dja_ 

dv„. 



dgp 
dv„ 



(A8c) 
(A8d) 
(A8e) 
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